Holographic Ricci dark energy: Interacting model and cosmological constraints 
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We extend the holographic Ricci dark energy model to include some direct, non- 
, gravitational interaction between dark energy and dark matter. We consider three phe- 

T-H ■ nomenological forms for the interaction term Q in the model, namely, Q is taken propor- 

\ tional to the Hubble expansion rate and the energy densities of dark sectors (taken to be 

Pde, Pm, and pde + Pm, respectively). We obtain a uniform analytical solution to the three 
interacting models. Furthermore, we constrain the models by using the latest observational 
data, including the 557 Union2 type la supernovae data, the cosmic microwave background 
^5 \ anisotropy data from the 7-yr WMAP, and the baryon acoustic oscillation data from the 

ff^ ■ SDSS. We show that in the interacting models of the holographic Ricci dark energy, a more 

OA \ reasonable value of f^mo will be obtained, and the observations favor a rather strong coupling 

between dark energy and dark matter. 

O: 

U- Our universe is undergoing an accelerated expansion, which was discovered via the observations 

r~| . of distant type la supernovae (SN) [1], about a decade ago. Combined analysis of various observa- 

tional data points to a cosmological model that contains dark energy, an exotic cosmic component 
Q ■ with negative pressure. What is more, dark energy is dominating the evolution of the current 

+i , universe: it occupies about 73% of the total energy of the universe. However, currently, the nature 

^ \ of dark energy is still enigmatic, and so the revelation of the nature of dark energy raises one of 

the biggest challenges for the modern fundamental science. 
. The cosmological constant A, first introduced by Einstein in 1917 [2], is an important candidate 

^ \ for dark energy, because it can provide a nice explanation for the accelerating universe and can 

fit the observational data well. Nevertheless, the cosmological constant is suffering from severe 

■ theoretical challenges [3] : one cannot explain why the theoretical value of A from the current 
CN \ framework of physics is greater than the observational value by many orders of magnitude, and 
CN ■ why the densities of dark energy and matter are in the same order just today while they evolve 

■ with rather different ways. Besides the cosmological constant, there is a fairly attractive idea that 
. the dark energy may not be a rigescent constant but a dynamically evolving component. In fact, 
I many dynamical models of dark energy have been proposed and studied in detail. Among the many 

• j ■ dynamical models of dark energy, the holographic model of dark energy [4-6] is very attractive 

. since it originates from the consideration of the holographic principle [7] of quantum gravity. It is 

\ widely believed that the physical nature of dark energy is in deep connection with the underlying 

quantum gravity theory. Thus, the theoretical and phenomenological studies on holographic dark 
energy may provide significant clues for the bottom-up exploration of a full quantum theory of 
gravitation. 

When considering gravity in a quantum field system, the conventional local quantum field theory 
will break down due to the too many degrees of freedom that would cause the formation of a black 
hole. According to the holographic principle, one may put an energy bound on the vacuum energy 
density, PvacL^ < Mp^L [8], where pvac is the vacuum energy density and Mpi is the reduced Planck 
mass. This bound implies that the total energy in a spatial region with size L should not exceed the 
mass of a black hole with the same size. The largest size compatible with this bound is the infrared 
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(IR) cutoff size of this effective field theory. Evidently, the holographic principle gives rise to a dark 
energy model basing on the effective quantum field theory with a UV/IR duality. From the UV/IR 
correspondence, the UV problem of dark energy can be converted into an IR problem. A given IR 
scale can saturate that bound, and so one can write the dark energy density as pde = 3c^MpiL~^ [4], 
where c is a phenomenological parameter (dimensionless) characterizing all of the uncertainties of 
the theory. Note that, hereafter, we use pde to denote the dark energy density. Now, the problem 
becomes how to choose an appropriate IR cutoff for the theory. Li [4] proposed to choose the event 
horizon of the universe as the IR cutoff of the theory. This choice generates a successful holographic 
model of dark energy, explaining both the fine-tuning problem and the cosmic coincidence problem 
at the same time, in some degree. Subsequently, other holographic models of dark energy basing 
on the different choices of IR cutoff were also proposed. For example, the choice of the conformal 
age of the universe leads to the agegraphic dark energy model [9] , and the choice of the Ricci scale 
of the universe gives rise to the holographic Ricci dark energy model [10]. 

The present paper focuses on the holographic Ricci dark energy (hereafter, abbreviated as RDE) 
model. In the RDE model, the IR length scale L is given by the average radius of the Ricci scalar 
curvature \TZ\ ' , so in this case the density of the holographic dark energy is p^c oc TZ. In a 
spatially flat universe, the Ricci scalar of the spacetime is given hy TZ = -6{H + 2H'^), where 
H = d/a is the Hubble expansion rate of the universe, and the dot denotes the derivative with 
respect to the cosmic time t. Therefore, the density of dark energy can be expressed as [10] 

Pde = 3aM^^{H + 2H^), (1) 

where a is a dimensionless parameter replacing of the Li model [4]. Note that, throughout the 
paper, we consider a flat universe owing to the fact that the flatness of the space is an important 
prediction of the inflationary cosmology and has been confirmed by observations, e.g., the current 
constraint is Qko ~ 10~^ [11]. So, the Friedmann equation is written as SMp^H^ = Pde + Pnu where 
Pm denotes the matter density. The RDE model has been studied extensively. In Ref. [12], the 
holographic meaning of RDE was revealed by investigating the causal connection scale -Rcc- In 
Ref. [13], the cosmological constraints on the RDE were studied, and the quintom feature was found. 
In Ref. [14], it was shown that the existence of the cosmic doomsday in the RDE model would ruin 
the theoretical foundation of the scenario, and a mend from the consideration of extra-dimension 
effects could erase the big-rip singularity and leads to a de Sitter finale for the holographic cosmos. 
For other relevant work, see, e.g., [15]. 

It is natural to consider the possible interaction between dark energy and dark matter in the 
RDE model. If dark energy interacts with cold dark matter, the continuity equations for them are 

Pdc + 3H{l + w)pdc = -Q, (2) 

Pm + 5Hpm = Q, (3) 

where w is the equation of state (EOS) parameter of dark energy, and Q is the interaction term. 
Note that, although p^ includes densities of cold dark matter and baryon matter, in this place we 
use pm to approximately describe dark matter density due to the fact that the density of baryon 
matter is much less than that of dark matter. In addition, owing to the lack of mechanism of 
microscopic origin of the interaction, one has to assume the forms of Q phenomenologically. For 
interacting dark energy model, several forms for Q have been put forward [16]. The most widely 
used form is Q oc Hp, where p denotes the energy density of the dark sectors, and usually it 
has three choices, namely, p = pdc, p = Pm, and p = pdc + Pm- Suwa and Nihei [17] considered 
the case of Q oc Hpdc for the interacting Ricci dark energy (IRDE) model. And, they placed the 
cosmological constraints on this model by using the SN, cosmic microwave background (CMB), and 
baryon acoustic oscillation (BAO) data. However, it should be pointed out that, in the work of 
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Suwa and Nihei [17], only one special case of interaction is considered. Moreover, the observational 
data used in the cosmological constraints are outdated today: 307 SN data from Union dataset, 
and CMB shift parameter from WMAP 5-year observation. In this paper, we will make some 
improvements. We will consider a more general interaction term, Q oc Hp, in the IRDE model, 
and place the cosmological constraints on the model by employing the latest observational data. 
Specifically, the interaction term is written as Q = 3bHp, where p denotes pde, Pm, and pde + Pm, 
respectively, and 6 is a dimensionless coupling constant. According to our convention, b > means 
that the energy transfer is from dark energy to cold dark matter. 

Combining with Eqs. (1) and (3), the Friedmann equation can be expressed as 



2 dx^ 
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where C = 1 — 7a/2 — 2>ah/2 and D = 3 — 6a — Qab for Q oc -ffpdci C = 1 — 7a/2 + 3a6/2 and 
D = Qa — Qah + 36 — 3 for Q oc H p^; C = 1 — 7a/2 and Z) = 3 — 6a — 36 for Q oc H{pde + Pm)', 
and X = Ina with a being the scale factor of the universe. The solution to Eq. (4) is obtained. 
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The initial conditions of Eq. (4) are 
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where OdeO = 1 ~ ^mo- The constants Aj- are given by 
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In what follows, we will constrain the IRDE model by using the latest observational data. 
We win use the 557 SN data from the Union2 dataset, the CMB data from the WMAP 7-year 
observation, and the BAO data from the SDSS. We will obtain the best-fitted parameters and 
likelihoods by using a Markov Chain Monte Carlo (MCMC) method. 

We use the data points of the 557 Union2 SN compiled in Ref. [18]. The theoretical distance 
modulus is defined as 



Pthizi) = 51ogio DL{zi) + po, 



(10) 



where z = l/a — lis the redshift, /io = 42.38 — 51ogio h with h the Hubble constant Hq in units of 
100 km s~^ Mpc~^, and the Hubble-free luminosity distance 



Dl{z) = {l + z) 



dz' 



(11) 
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where E{z) = H{z)/Ho, and 6 denotes the model parameters. Correspondingly, the function 
for the 557 Union2 SN data is given by 



>^sn(^) = 2^ , (12) 



i=l 

where a is the corresponding la error of distance modulus for each supernova. The parameter //q 
is a nuisance parameter and one can expand Eq. (12) as 

xhie) = A{9) - 2f,oBie) + f^lC (13) 

where A{6), B{6) and C are defined in Ref. [19]. Evidently, Eq. (13) has a minimum for /^o = B /C 
at 

xUe) = A{e) - (14) 

Since Xsn min ~ XsN min' iiistead minimizing xIn ^^^^ minimize Xsn which is independent of 
the nuisance parameter fiQ. 

Next, we consider the cosmological observational data from WMAP and SDSS. For the WMAP 
data, we use the CMB shift parameter R; for the SDSS data, we use the parameter A of the BAO 
measurement. It is widely believed that both R and A are nearly model-independent and contain 
essential information of the full WMAP CMB and SDSS BAO data [20]. The shift parameter R is 
given by [20, 21] 



R^^ll! I ^> (15) 



dz 

where the redshift of recombination z^ = 1091.3, from the WMAP 7- year data [11]. The shift 
parameter R relates the angular diameter distance to the last scattering surface, the comoving 
size of the sound horizon at and the angular scale of the first acoustic peak in the CMB power 
spectrum of temperature fluctuations [20, 21]. The value of R is 1.725 it 0.018, from the WMAP 
7- year data [11]. On the other hand, the distance parameter A from the measurement of the BAO 
peak in the distribution of SDSS luminous red galaxies [22] is given by 



A ^ n]l^E{z^)-'/-' 



1 dz 



Zh Jo E{z) 



2/3 

(16) 



where z^ = 0.35. In Ref. [23], the value of A has been determined to be 0.469 (7is/0.98)-°-^^ ±0.017. 
Here, the scalar spectral index ng is taken to be 0.963, from the WMAP 7- year data [11]. So the 
total is given by 

= xiN + XcMB + xlAO' (17) 

where Xsn is given by (14), Xcmb = (-^ - Rohsf/cr% and x|ao = (^ " ^obs) Vo"!- The best-fitted 
model parameters are determined by minimizing the total y^. 

Now we fit the RDE and IRDE models to the observational data. We use the MCMC method 
and finally we obtain the best-fit parameters and the corresponding Xmin- The best-fit, la and 2a 
values of the parameters with Xmin ^'^^'^ models are all presented in Table I. For convenience, 

in the table we abbreviate the three interacting models: IRDEl corresponds to Q = 36-ffpde) 
IRDE2 corresponds to Q = 3bHp^, and IRDE3 corresponds to Q = 3bH{pdc + Pm)- 
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TABLE I: The fit results of tlic ACDM, RDE and IRDE models. Here, IRDEl model corresponds to 
Q = 36-ffpdo, IRDE2 corresponds to Q = 3bHpn^, and IRDE3 model corresponds to Q = 3bH{pdc + Pm)- 



Model 



n 



mO 



a 



AAIC ABIC 



-0.013-0.026 
n QOQ+0. 022+0. 037 
U.OZO_Q 021-0.035 
n 970+0.034+0.051 
U.Z/0_o. 032-0. 048 
n OfiQ+0. 034+0. 052 
U.ZO»_Q 032-0.049 
hO. 034+0. 052 
-0.032-0.048 



ACDM 
RDE 
IRDEl 
IRDE2 
IRDE3 



0.270^ 



0.270: 



N/A 
N/A 



0.052 
0.032 
0.020 



+0.011+0.017 
0.015-0.025 
+0.013+0.021 
0.013-0.019 
0.006+0.010 
0.007-0.011 



N/A 



n Qr;f;+0. 021+0. 035 
U.OOU_o. 021-0. 035 
n /I/I9+0. 045+0. 069 
"•^^^-0. 044-0. 066 
n /IQQ+0.043+0.066 
U.'iOO_o 042-0.062 
n /1 07+0. 045+0. 068 
U.^J/_0 043-0.064 



542.919 
565.683 
542.698 
543.038 
542.875 





24.764 
3.779 
4.119 
3.956 





29.090 
12.431 
12.771 
12.608 




0.28 0.30 0.32 0.34 0.36 

n „ 

mO 

FIG. 1: The probability contours at Icr and 2a confidence levels in the f^mo ~ " plane for the RDE model. 




FIG. 2: The probability contours at la and 2a confidence levels in the ft^Q — a and b — a planes for the 
IRDEl model (corresponding to Q = SbHp^e)- 
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FIG. 4: The probability contours at la and 2cr confidence levels in the ft^o ^ ol and & — a planes for the 
IRDE3 model (corresponding to Q = 3bH{pdc + Pm))- 



For the RDE model, we obtain the best-fit values of the parameters: OmO = 0.323 and a = 0.356, 
and the corresponding minimal is Xmin ~ 565.683. Obviously, the fit value of f^mo in the RDE 
model is rather big, evidently greater than that of the Lambda Cold Dark Matter (ACDM) model, 
0.270. Note that for making a comparison with the ACDM model, we also fit the ACDM model 
and list the results also in Table I. Now, let us see the fit results of the IRDE models. We find 
that, for all the three IRDE models, the best-fit values obtained for r^mo and a are approximately 
equal: OmO = 0.27 and a = 0.44 (for IRDE2, a = 0.43). So, it is clearly seen that the fit value of 
in the IRDE models is more reasonable, consistent with that of the ACDM model. Comparing 
the values of a, it is found that, when there is an interaction between dark energy and cold dark 
matter, the value of a tends to become bigger. From the above analysis on the fitting results, we 
find that the parameters fimo a-nd a should be in anti-correlation. For the coupling constant b in 
the IRDE models, we see clearly that the values of b are explicitly different for the three IRDE 
models. The value of b in the IRDEl model is the biggest, b = 0.052; the value of b in the IRDE3 
model is the smallest: b = 0.020. Since in the IRDE models there is an additional parameter, b, 
the values of Xmin should be substantially less than that of the RDE model. This is the case: the 
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three IRDE models have the similar values of Xmin' around 543, much smaller than that of the 
RDE model. Among the three IRDE model, the IRDEl model has the smallest Xmin' 542.698. 

However, it is unwise to use the statistic solely to compare competing models since the 
number of parameters is different for the models. In general, a model with more parameters tends 
to give a lower Xmin' may employ the information criteria (IC) to assess the worth of 

models. These statistics favor models that give a good fit with fewer parameters. Thus, we shall 
use the AIC (Akaike information criterion) and BIC (Bayesian information criterion) as model 
selection criteria, defined, respectively, as 

AIC = Xmin + 2A:, BIC = xLn + klnN, (18) 

where k is the number of parameters, and is the data points used in the fit. Note that the 
absolute value of the criterion is not of interest, only the relative value between different models, 
AAIC or ABIC, is useful. For the details about the AIC and BIC, especially their applications 
in a cosmological context, see, e.g., Refs. [24-28]. We give the calculation results of AIC and BIC 
for the models in Table I. Since the ACDM model is an important reference model for the studies 
of dark energy models, we also include its results in the table. From Table I, we see that, though 
the value of Xmin ACDM model is greater than that of the IRDE models, the ACDM model 

gives the lowest values of AIC and BIC. This result is consistent with the previous studies [26-28]. 
Actually, the ACDM model is the best model fitting the current data among all of the existing dark 
energy models. Therefore, the ACDM model is chosen to be a fiducial model when we present the 
values of AAIC and ABIC. Obviously, the RDE model is strongly unsupported by the data, since 
its ABIC is fairly high; see also Ref. [27] for comparison of various dark energy models including 
the RDE model. Comparing with the RDE model, the IRDE models are much better, though they 
are still worse than the ACDM model. 

Figures 1—4 show the probability contours at 68.3% and 95.4% confidence levels in — " 
and b — a parameter planes for the RDE and IRDE models. We find that the r^mo — planes of 
the IRDE models show a stronger degeneracy than that of the RDE model. In addition, for all 
the three IRDE models, we see that J^mO and a are anti-correlated, and b and a are in positive 
correlation. 

It should also be noted that, in the work of Suwa and Nihei [17], the contour plots in the 
parameter planes are incomplete: when plotting the a — r^deo contours, 7 = 36 is fixed to be 0.15; 
when plotting the 7 — 0: contours, f^deo is fixed to be 0.73; see Figs. 2 and 3 in Ref. [17]. In our 
work, the probe of the parameter-space is complete, i.e., when plotting the contours in parameter 
plane, other nuisance parameters are marginalized. In addition, it should be mentioned that the 
IRDEl case was also investigated in Ref. [28]. In Ref. [28], the authors analyzed in detail the 
various holographic dark energy models, with the IR cutoff chosen to be the Hubble scale, the 
future event horizon, and the Ricci scale, respectively, with an interaction between dark energy 
and dark matter. For the IRDEl model, they also constrain the model by using the current data, 
but they assume wq = —1 in their analysis. The result of the coupling constant they obtained is 
b = 0.05 lb 0.01, almost the same as ours. Therefore, we find that our study confirms the results of 
Ref. [28]. For other relevant work, see also, e.g., Ref. [29]. 

In summary, we have extended the RDE model to include the possible interaction between dark 
energy and cold dark matter in this paper. Though in Refs. [17, 28] a specific form of the interaction 
in the RDE model has been investigated, our study makes some significant improvements. In 
this work, we have considered three phenomenological forms for Q in the RDE model, namely, 
Q = 3bHp, where p can be taken to be pdci Pm 

and pdc + Pm, respectively. We have solved the 
IRDE models, and got a uniform analytical solution for the three IRDE models. Furthermore, we 
constrained the RDE and IRDE models by using the latest observational data, including the 557 
Union2 SN data, the CMB WMAP 7-yr data, and the BAG SDSS data. Our fit results show that 
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when the interaction between dark energy and dark matter is considered in the RDE model, a more 
reasonable value of OmO will be obtained, i.e., flmo = 0.27, rather than the value of fimo = 0.32 in 
the RDE model without interaction. Moreover, it has been shown that a rather strong coupling 
between dark energy and cold dark matter in the RDE model is favored by the observations: the 
best-fit value of b is around 0(10"^), in particular, for the IRDEl model, b = 0.052. Our work 
makes the study of the interacting model of the holographic Ricci dark energy more complete. 
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